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ABSTRACT 

The backpressure routing policy is known to be a throughput 
optimal policy that supports any feasible traffic demand in 
data networks, but may have poor delay performance when 
packets traverse loops in the network. In this paper, we 
study loop-free backpressure routing policies that forward 
packets along directed acyclic graphs (DAGs) to avoid the 
looping problem. These policies use link reversal algorithms 
to improve the DAGs in order to support any achievable 
traffic demand. 

For a network with a single commodity, we show that a 
DAG that supports a given traffic demand can be found after 
a finite number of iterations of the link-reversal process. We 
use this to develop a joint link-reversal and backpressure 
routing policy, called the loop free backpressure (LFBP) 
algorithm. This algorithm forwards packets on the DAG, 
while the DAG is dynamically updated based on the growth 
of the queue backlogs. We show by simulations that such 
a DAG-based policy improves the delay over the classical 
backpressure routing policy. We also propose a multicom¬ 
modity version of the LFBP algorithm, and via simulation 
we show that its delay performance is better than that of 
backpressure. 

1. INTRODUCTION 

Throughput and delay are the two major metrics used to 
evaluate the performance of communication networks. For 
networks that exhibit high variability, such as mobile ad hoc 
networks, the dynamic backpressure routing policy [1] is a 
highly desirable solution, known to maximize throughput in 
a wide range of settings. However, the delay performance 
of backpressure is poor [2]. The high delay is attributed to 
a property of backpressure that allows the packets to loop 
within the network instead of moving towards the destina¬ 
tion. In this paper we improve the delay performance of 
backpressure routing by constraining the data routing along 
loop free paths. 

To eliminate loops in the network, we assign directions to 
the links such that the network becomes a directed acyclic 
graph (DAG). Initially, we generate an arbitrary DAG and 
use backpressure routing over it. If the initial DAG has max- 
flow smaller than the traffic demand, parts of the network 


become overloaded. By reversing the direction of the links 
that point from non-overloaded to overloaded nodes a new 
DAG with a lower overload is obtained. Iterating over this 
process, our distributed algorithm gradually converges to a 
DAG that supports any traffic demand feasible in the net¬ 
work. Hence the loop-free property is achieved without the 
loss of throughput. 

Prior work identifies looping as a main cause for high de¬ 
lays in backpressure routing and proposes delay-aware back¬ 
pressure techniques. Backpressure enhanced with hop count 
bias is first proposed in [3] to drive packets through paths 
with smallest hop counts when the load is low. An alter¬ 
native backpressure modification that utilizes shortest path 
information is proposed in [8]. A different line of works pro¬ 
poses to learn the network topology using backpressure and 
then use this information to enhance routing decisions. In 
[7] backpressure is constrained to a subgraph which is dis¬ 
covered by running unconstrained backpressure for a time 
period and computing the average number of packets routed 
over each link. Learning is effectively used in scheduling [9] 
and utility optimization [13] for wireless networks. In our 
work we aim to eliminate loops by restricting backpressure 
to a DAG, while we dynamically improve the DAG by re¬ 
versing links. 

The link-reversal algorithms were introduced in [4] as a 
means to maintain connectivity in networks with volatile 
links. These distributed algorithms react to any topological 
changes to obtain a DAG such that each node has a loop-free 
path to the destination. In [5], one of the link-reversal algo¬ 
rithms was used to design a routing protocol (called TORA) 
for multihop wireless networks. Although these algorithms 
provide loop free paths and guarantee connectivity from the 
nodes to the destination, they do not maximize through¬ 
put. Thus, the main goal of this paper is to create a new 
link-reversal algorithm and combine it with the backpressure 
algorithm to construct a distributed throughput optimal al¬ 
gorithm with improved delay performance. 

The main contributions of this paper are as follows: 

• For a DAG with a single commodity, we study the lex¬ 
icographic optimization of the queue growth rate. We 
develop a novel link-reversal algorithm that reverses 
link direction based on overload conditions to form a 



new DAG with lexicographically smaller queue growth 
rates. 

• We show that the queue backlog information from back¬ 
pressure routing can be used to perform the link re¬ 
versals. We develop loop free backpressure (LFBP) al¬ 
gorithm, a distributed routing scheme that eliminates 
loops and retains the throughput optimality property. 

• Our simulation results of LFBP show a signihcant de¬ 
lay improvement over backpressure in static and dy¬ 
namic networks. 

• We extend the LFBP algorithm to networks with mul¬ 
tiple commodities, and provide a simulation result to 
show its delay improvement over backpressure. 

2. SYSTEM MODEL AND DEFINITIONS 

2.1 Network model 

We consider the problem of routing single-commodity data 
packets in a network. The network is represented by a graph 
G = {N, E), where N is the set of nodes and E is the set of 
undirected links {i,j} with capacity dj. Packets arrive at 
the source node s at rate A and are destined for a receiver 
node d. Let denote the maximum flow from node s 

to d in the network G. The quantity is the maximally 
achievable throughput at the destination node d. 

To avoid unnecessary routing loops, we restrict forwarding 
along a directed acyclic graph (DAG) embedded in the graph 
G. An optimal DAG exists to support the max-flow 
and can be found by: (i) computing a feasible flow allocation 
(fij) that yields the max-flow in G (e.g. using [11]); (ii) 
trimming any positive flow on directed cycles; (iii) defining 
an embedded DAG by assigning a direction for each link 
{i,j} according to the direction of the flow fij on that link. 
Since backpressure achieves the max-flow of a constrained 
graph [14], performing backpressure routing over the optimal 
DAG supports A. 

This centralized approach is unsuitable for mobile ad- 
hoc networks, which are based on wireless links with time- 
varying capacities and may undergo frequent topology changes. 
In such situations, the optimal embedded DAG also changes 
with time, which requires constantly repeating the above 
offline process. Instead, it is possible to use a distributed 
adaptive mechanism that reverses the direction of links un¬ 
til a DAG that supports the current trafhc demand is found. 

In this paper we propose an algorithm that reacts to the traf¬ 
fic conditions by changing the directions of some links. To 
understand the properties of the link-reversing operations, 
we first study the fluid level behavior of a network under 
overload conditions. 

2.2 Flow equations 

Consider an embedded DAG Dk = {Nk,Ek) in the net¬ 
work graph G, where A^fc = is the set of network nodes 
and Ek is the set of directed links.^ For each link {i,j} € E, 
either {i,j) or {j,i) belongs to Ek (but not both). Each 
directed link {i,j) has the capacity of the undirected coun¬ 
terpart {i,j}, which is aj. Let be the maximum flow 

of the DAG Dk from the source node s to the destination 

^The notation Dk of an embedded DAG is useful in the 
paper; it will denote the DAG that is formed after the fcth 
iteration of the link-reversal algorithm. 


node d. Any embedded DAG has smaller or equal max-flow 
with respect to G, /(”“ < 

For two disjoint subsets A and B of nodes in Dk, we define 
capj,(A, B) as the total capacity of the directed links going 
from A to B, i.e., 

capk{A,B)= (1) 

(i,3)eEk:ieA, jSB 

A cut is a partition of nodes {A, A’^) such that s £ A 
and d £ A'^. A cut {Ak,A%) is a min-cut if it minimizes 
the expression capj,(Afe, A^) over all cuts. By the max-flow 
min-cut theorem = capj,(Afc, A^), where {Ak,A%) is 

the min-cut of the DAG Dk- We remark that a cut in a 
DAG is also a cut in G or another DAG. However, the value 
of capj,(.,.) depends on the graph considered (see summation 
in (1)). 

We consider the network as a time-slotted system, where 
slot t refers to the time interval [t, t -|- 1), t £ {0,1,2,...}. 
Each network node n maintains a queue Qnit), where Qn{t) 
also denotes the queue backlog at time t. We have Qd(t) = 0 
for all t since the destination node d does not buffer packets. 
Let A{t) be the number of exogenous packets arriving at the 
source node s in slot t. Under a routing policy that forwards 
packets over the directed links defined by the DAG Dk, let 
Eij (t) be the number of packets that are transmitted over the 
directed link {i,j) £ Ek in slot t; the link capacity constraint 
states that Fij{t) < dj for all t. The queues Qn{t),n ^ d, 
are updated over slots according to 

(3n(t)=gn(t-l) + l[„=s]4l(t) 

+ E E ( 2 ) 


where ![.] is an indicator function. 

To study the overload behavior of the system we define 
the queue overload (i.e., growth) rate at node n as 

,. Qnii) 

q„ = hm —-—. (3) 

t—¥00 t 

Additionally, dehne the exogenous packet arrival rate A and 
the flow fij over a directed link {i,j) as 


A = lim 


E^M’ 


fij = lim -E-^b('r), 

t—foo t ^^ 

T—0 


where the above limits are assumed to exist almost surely 
(see [6] for details). Using the recursion (2), taking time 
averages and letting t —>■ oo, we have the fluid-level equation: 


Qn — l[n—s]A-t- 'y ] fin ^ ^ fnj, Vu £ A^\{d} 

i:(i,n)eEt, j-{n,j)eEk 

(4) 

0 < fij < dj, y{i,j) e Ek- (5) 

Equations (4) and (5) are the flow conservation and link ca¬ 
pacity constraints, respectively. A network node n is said to 
be overloaded if its queue growth rate Qn is positive, which 
implies that Q„{t) —>■ cxa as t —>■ oo (see (3) and [10]). Sum¬ 
ming (4) over n £ N yields 


E qn = \- E 

n^N\{d} 



where X]i (i d)6Bfc denotes the throughput received at the 
destination d. Therefore, equation (6) states that the re¬ 
ceived throughput is equal to the exogenous arrival rate A 
less the sum of queue growth rates network. 

2.3 Properties of queue overload vector 

If the traffic arrival rate A is strictly larger than the max¬ 
imum flow of the DAG Dk, then some network nodes 
will be overloaded. It is because, from (6), we have 

= /w>A-/r">o, (7) 

nGN i:(i,d)£E^ 

which implies that Qn > 0 for some node n £ N. Let 
q = {qn)neN be the queue overload vector. A queue overload 
vector q is feasible in the DAG Dk if there exist overload 
rates {qn)neN and flow variables {fij)(i,j)£Ek that satisfy (4) 
and (5). Let Qk be the set of all feasible queue overload vec¬ 
tors in Dk- We are interested in the lexicographically small¬ 
est queue overflow vector in set Qk- Formally, given a vector 
u = (ui, - - - ,un), let Mi be the ith maximal component of 
u- We say that a vector u is lexicographically smaller than 
a vector v, denoted by u <iex u, if < ui or Ui = Vi for 
all i = 1,..., {j — 1) and Uj < Vj for some j = 2, - - -, N- If 
Ui = Vi for all i, then the two vectors are lexicographically 
equal, represented by u =iex v-^ The above-defined vector 
comparison induces a total order on the set Qk, and hence 
the existence of a lexicographically smallest vector is always 
guaranteed [12]. 

Lemma 1 ([6]). Let q™” be the lexicographically small¬ 

est vector in the queue overload region Qk of the DAG Dk- 
We have the following properties: 

1- The vector q™" exists and is unique in the set Qk- 

2- The vector q™" minimizes the sum of queue overload 
rates, i.e., it is a solution to the optimization problem: 

minimize subject to q £ Qk 

nSN 

(direct consequence of Theorem 1 in [6])- Due to (6), 
the corresponding throughput is maximized- 

3- A feasible flow allocation vector ifij)(i,j)eEk induces 
q™" if and only if over each link {i,j) £ Ek the fol¬ 
lowing holds: 

if qi < qj, then flj = 0, (8) 

if T > Qj, then flj = dj- (9) 

In general, there are many flow allocations that yield the 
maximum throughput. Focusing on those that additionally 
induce q™‘" has two advantages. First, as shown next, these 
allocations lead to link-reversal operations that improve the 
max-flow of the DAG Dk- Second, the backpressure algo¬ 
rithm can be used to preform the same reversals and improve 
the max-flow; we will use this observation in Section 4 to 
combine link-reversal algorithms with backpressure routing. 

^As an example, the two vectors u = (3, 2,1, 2,1) and v — 
(1,2,3,2,2) satisfy u <iex v because ui = ui = 3, U 2 = 
V2 = U3 = V3 = 2, and U4 = 1 < U4 = 2. 


3. LINK-REVERSAL ALGORITHMS 

The link-reversal algorithms given in [4] were designed to 
maintain a path from each node in the network to the des¬ 
tination. One algorithm relevant to this paper is the full re¬ 
versal method- This algorithm is triggered when some nodes 
n ^ d lose all of their outgoing links. At every iteration of 
the algorithm, nodes n, that have no outgoing link, reverse 
the direction of all their incoming links. This process is re¬ 
peated until all the nodes other than the destination have at 
least one outgoing link. When the process stops these nodes 
are guaranteed to have a path to the destination. The ex¬ 
ample in Figure 1, taken from [4], illustrates this algorithm 
at work. 



(d) (e) 

Figure 1: Illustration of the full reversal method of 
[4] when the dashed link in Figure 1(a) is lost. At 
every iteration, the algorithm reverses all the links 
incident to the nodes with no outgoing link. 

Although the full reversal algorithm guarantees connectiv¬ 
ity, the resulting throughput may be significantly lower than 
the maximum possible. Hence, in this paper we shift the 
focus from connectivity to maximum throughput. Specifi¬ 
cally, we propose a novel link-reversal algorithm that pro¬ 
duces a DAG which supports the traffic demand A, assum¬ 
ing A < We do this by quickly constructing an initial 

DAG and improving upon it in multiple iterations. 

3.1 Initial DAG 

We assume that each node in the network has a unique 
ID. These IDs give a topological ordering to the nodes. So, 
the initial DAG can be created simply by directing each link 
to go from the node with the lower ID to the node with the 
higher ID. If the unique IDs are not available, the initial 
DAG can be created by using a strategy such as the one 
given in [5]. 

3.2 Overload detection 

Given a DAG Dk, k — 0,1,2, - - -, we suppose that there is 
a routing policy tt that yields the lexicographically minimal 
queue overload vector q™‘“. ^ Then we use the vector q™‘" 
to detect node overload and decide whether a link should be 
reversed. 

^Such a policy tt can simply solve an optimization problem 
offline to compute the required flow allocation. In Section 4, 
we develop a distributed algorithm using backpressure that 
does not require the computation of the lexicographically 
optimal overload vector. We use this vector only to prove 
the properties of our link-reversal algorithm. 






If the data arrival rate A is less than or equal to the max¬ 
imum flow of the DAG Dk, then there exists a flow 

allocation (fij) that supports the traffic demand and yields 
zero queue overload rates q„ = 0 sd. all nodes n £ N. By the 
second property of Lemma 1 and nonnegativity of the over¬ 
load vector, the queue overload vector is zero. Thus, 
the throughput under policy tt is A according to (6), and the 
current DAG Dk supports A; no link-reversal operations are 
needed. 

On the other hand, if the arrival rate A is strictly larger 
than the maximum flow by the second property in 

Lemma 1 the maximum throughput is and the queue 

overload vector = (9™n)neiv is nonzero because we 

have from (7) that 

•Jfe.n > A - /fe >0. 

nSiV 

We may therefore detect the event “DAG Dk supports A” by 
testing whether the overload vector q™'" is zero or non-zero. 

The next lemma shows that if DAG Dk does not support 
A then it contains at least one under-utilized link (our link- 
reversal algorithm will reverse the direction of such links to 
improve network throughput). 

Lemma 2. Suppose that the traffic demand A satisfies 

< A < /”“. 

where is the max-flow of the DAG Dk and js the 
max-flow of the undirected network G. Then there exists a 
link {i,j) € Ek such that = 0 and q™f > 0 . 

Proof of Lemma 2. Let Ak be the set of overloaded 
nodes under a flow allocation that induces the lexicograph¬ 
ically minimal overload vector in the DAG Dk', the set 
Ak is nonempty due to A > and (7). It follows that the 
partition {Ak,Al) is a min-cut of Dk (see Lemma 7 in the 
Appendix)."* By the max-flow min-cut theorem, the capac¬ 
ity of the min-cut {Ak,Ak) in Dk satisfies capki^k, Af) = 

ymax ^ ymax 

The proof is by contradiction. Let us assume that there is 
no directed link that goes from the set Af to Ak in the DAG 
Dk- It follows that cap k{Ak, Af) is the sum of capacities of 
all undirected links between the sets Ak and Af, i.e., 

capk{Ak,Al)= Cij, 

ieAk,j^Ak 

which is equal to the value of the cut {Ak, Af) in graph G. 
Since the value of any cut is larger or equal to the min-cut, 
applying the max-flow min-cut theorem on G we have 

< 2 ^ Cij = cap k{Ak,Ak) = fk , 

iCiAk , j^Ak 

which contradicts the assumption that < A < □ 

3.3 Link reversal 

Lemma 2 shows that if the DAG Dk has insufficient capac¬ 
ity to support the traffic demand A < then there exists 

a directed link from an underloaded node i to an overloaded 
one j under the lexicographically minimum overflow vector 
q™'". Because of property (8), we may infer that this link 

^The set Af contains the destination node d and is 
nonempty. 


is not utilized. Next we show that reversing the direction of 
this link provides a strictly improved DAG. 

We consider the link-reversal algorithm (Algorithm 1) that 
reverses all such links that satisfy the property in Lemma 2. 
This reversal yields a new directed graph Dk-ei ~ {N, Ek+i)- 


Algorithm 1 Link-Reversal Algorithm 
1: for all {i,j) £ Ek do 
2: if q^ffl = 0 and q^^ > 0 then 

3: {j, %) £ Ek-\-\ 

4: else 

5: {i,j)£Ek+\ 

6: end if 

7: end for 


Lemma 3. The directed graph Dk+i is acyclic. 

Proof of Lemma 3. Recall that Ak is the set of over¬ 
loaded nodes in the DAG Dk under the lexicographically 
minimum queue overload vector q™‘“. Let Lk C E he the 
set of undirected links between Ak and A%. Algorithm 1 
changes the link direction in a subset of Lk- More precisely, 
it enforces the direction of all links in Lk to go from Ak to 

We complete the proof by construction in two steps. First, 
we remove all links in Lk from the DAG Dk, resulting in two 
disconnected subgraphs that are DAGs themselves. Second, 
consider that we add a link in Lk back to the network with 
the direction going from Ak to Af. This link addition does 
not create a cycle because there is no path from Af. to Ak, 
and the resulting graph remains to be a DAG. We can add 
the other links in Lk one-by-one back to the graph with the 
direction from Ak to A%; similarly, these link additions do 
not create cycles. The final directed graph is Dk+i, and it 
is a DAG. See Fig. 2 for an illustration. □ 

The next lemma shows that the new DAG Dk+i sup¬ 
ports a lexicographically smaller optimal overload vector 
(and therefore potentially better throughput) than the DAG 

Dk- 

Lemma 4. Let Dk be a DAG with the maximum flow 
< A < The DAG Dk-ei, obtained by perform¬ 

ing Algorithm 1 over Dk, has the lexicographically minimum 
queue overload vector satisfying qffff <iex q™”- 

Proof of Lemma 4. Gonsider a link (a, 6) e Ek such 
that = 0 and q^ff > 0; this link exists by Lemma 2. 
From the property (8), any feasible flow allocation (fij) that 
yields the lexicographically minimum overload vector q™‘" 
must have fab = 0 over link (a, fe). The link-reversal algo¬ 
rithm reverses the link (a, 6) so that {b, a) £ Ek+i in the 
DAG Dk+i- Consider the following feasible flow allocation 
{fij) on the DAG Dk+p- 

(e if (i, j) = (fe,a) 

fij = lo = fji if {i,j) f {b,a) but {j,i) is reversed 
[fij if {i,j) is not reversed 

where e < q™)," is ^ sufficiently small value. In other words, 
the flow allocation {fij) is formed by reversing links and 
keeping the previous flow allocation {fij) except that we 





(a) The DAG Dk with (b) Two disconnected 
Aj; = {s, 2, 3, 5}. DAGs formed by re¬ 
moving all links be¬ 
tween Ak and A%. 



(c) The DAG Dk+i 
formed by adding all 
links in Lk back to the 
graph with the direc¬ 
tion going from Ak to 


Figure 2: Illustration for the proof of Lemma 3. 


0 .- = 

(a) Link {a, b) before the (b) Link (6, a) after the 
link reversal. link reversal. 


Figure 3: A link {a, 6} in the network in Fig. 2 before 
and after link reversal. Before the reversal, the flow 
fab is zero on (a, 6). After the reversal, an e flow can 
be sent over {b, a) so that (g^, 5t,) <iex while 

the rest of the flow allocation remains the same. 


forward an e-amount of overload traffic from node 6 to a. 
Let q = {qn)neN be the resulting queue overload vector. 
We have 


be distinct. Since there are a finite number of unique DAGs 
in the network, the link-reversal algorithm will find a DAG 
Dk* that has the lexicographically minimal overload vector 
“ = 0 and the maximum flow > A in a finite number 
of iterations; this DAG Dk» exists because the undirected 
graph G has the maximum flow > A. □ 

3.4 Arrivals outside stability region 

We show that even when A > the link reversal al¬ 

gorithm will stop reversing the links in a finite number of 
iterations, and it will obtain the DAG that supports the 
maximum throughput We begin by examining the 

termination condition of our algorithm and show that if the 
algorithm stops at iteration k, then the DAG Dk supports 
the max-flow of the network. 

Lemma 5. Consider the situation when A > If 

there is no link (i,j) such that = 0 and > 0, then 

jmax ^ jmax ^ ^max 

reverse at iteration k, and > 0, then the throughput of 
Dk is equal to 

Proof. Let Ak be the set of overloaded nodes under 
a flow allocation that induces the lexicographically mini¬ 
mal overload vector in the DAG Dk- We know that 
{Ak,Ak) is a min-cut of the network from Lemma 7 (in the 
appendix), so 

capk{Ak,Al) = fr^. 

Suppose the link reversal algorithm stops after iteration 
k, i.e. at iteration k there are no links to reverse. In this 
situation, there is no link (i, j) such that = 0 and > 
0, so by property (9), all the links between Ak and A| go 
from Ak to A^. The capacity of the cut {Ak, Af) is given by 

capk{Ak,Al)= Cij. 

ieAkJeAi 

This is equal to the capacity of the cut (A*,, A|) in the undi¬ 
rected network G. So f™^^ < capk{Ak, A%) = Be¬ 
cause /max cannot be greater than By 

assumption A > so A > □ 


min , mm mm r» i 

qb = qk,b - e < Qk^b , qa = e> qk,a = 0, and 
= qT'n, n ^ {a,b}. 

Therefore, q <iex <7™'" (see Fig. 3 for an illustration). Let 
(fk+i be the lexicographically minimal overload vector in 
Dk+i- It follows that q™/" <iex Q <iex q™'“, completing 
the proof. □ 

Theorem 1. Suppose the traffic demand is feasible in G, 
i.e., A < f™^^, and the routing policy induces the overload 
vector q™‘" at every iteration k. Then, the link-reversal al¬ 
gorithm will find a DAG whose maximum flow supports A in 
a finite number of iterations. 

Proof of Theorem 1. The link-reversal algorithm cre¬ 
ates a sequence of DAGs {Do, Di, D2, ...} in which a strict 
improvement in the lexicographically minimal overload vec¬ 
tor is made after each iteration, i.e., 

min .. min .. min .. 

Qo >1 ex Qi ^lex Q2 ^lex • • • . 

The lexicographically minimal overload vector is unique in 
a DAG by Lemma 1, the DAGs {Do, Di, D 2 ,...} must all 


When A > this lemma shows that the link reversal 

algorithm stops only when the DAG achieves the maximum 
throughput of the network. Hence, if the DAG doesn’t sup¬ 
port the maximum throughput, then there exists a link that 
can be reversed. After each reversal, Lemma 3 holds, so 
the directed graph obtained after the reversal is acyclic. We 
can modify Lemma 4 to show that every reversal produces 
a DAG that supports an improved lexicographically optimal 
overload vector. We can combine these results to prove the 
following theorem. 

Theorem 2. Suppose the traffic demand is not feasible 
in G, i.e., A > and the routing policy induces the 

overload vector q™‘" at every iteration k. Then, the link- 
reversal algorithm will find a DAG whose maximum flow 
supports in a finite number of iterations. 

4. DISTRIBUTED DYNAMIC ALGORITHM 

In the previous sections we developed a link reversal algo¬ 
rithm based on the assumption that we had a routing policy 
that lexicographically minimized the overload vector q™'". 






The algorithm reversed all the links that went from the set 
of all the non-overloaded nodes A'^ to the set of overloaded 
nodes Ak- We showed that repeating this process for some 
iterations results in a DAG that supports the arrival rate A. 

The goal of this paper is to develop a link reversal al¬ 
gorithm based on backpressure. To achieve this goal, we 
develop a threshold based algorithm that identifies the cut 
(Afe, Afc) using the queue backlog information of backpres¬ 
sure. We can use this cut to perform the same link rever¬ 
sals as in Section 3 without computing the lexicographically 
minimum overload vector. Hence, this algorithm generates 
the same sequence of DAGs as the link reversal algorithm 
described in the previous section. So, all the previous theo¬ 
rems hold, and the new algorithm will obtain the DAG that 
supports the arrival rate A (when possible). We will call this 
algorithm the loop free backpressure (LFBP) algorithm. 

We begin by creating an initial DAG Dq using the method 
presented in Section 3.1. Then, we use the backpressure al¬ 
gorithm to route the packets from the source to the desti¬ 
nation over Do- Let Q„{t) be the queue length at node n 
in slot t. The backpressure algorithm can be written as in 
Algorithm 2. It simply sends packets on a link {i,j) if node 
i has more packets than j. 


Algorithm 2 Backpressure algorithm (BP) 

1: for all {i,j) £ Ek do 
2: if Qi{t) > Qj{t) then 

3: Transmit mm{cij,Qi{t)} packets from i to j 

4: end if 

5: end for 


Since backpressure is throughput optimal [1], if the arrival 
rate is less than then all queues are stable. If the 

arrival rate is larger than the system is unstable and 

the queue length grows at some nodes. In this case, the 
next lemma shows that if we were using a routing policy 
that produced the optimal overload vector q™‘", the set of 
all the overloaded nodes Ak and the non-overloaded nodes 
A% form the smallest min-cut of the DAG Dk- 

Definition 1. We define the smallest min-cut {X*, X*^) 
in the DAG Dk as the min-cut with the smallest number of 
nodes in the source side of the cut, i.e., {X*,X*‘^) solves 

minimize: \X\ 

subject to: (X, A“) is a min-cut of Dk- 

Lemma 6. Let Ak the set of overloaded nodes under a 
flow allocation (fij) that induces the lexicographically min¬ 
imum overload vector in the DAG Dk- If \Ak\ > 0, then 
{Ak,A%) is the unique smallest min-cut in Dk- 

Proof of Lemma 6. The proof is in Appendix B. □ 

Essentially, at every iteration, the link reversal algorithm 
of Section 3 discovers the smallest min-cut {Ak,A%) of the 
DAG Dk and reverses the links that go from Aj. to Ak- 
Now the following theorem shows that the backpressure al¬ 
gorithm can be augmented with some thresholds to identify 
the smallest min-cut. 

Theorem 3. Assume that [Ak,A%) is the smallest min- 
cut for DAG Dk with a cut capacity of = cap{Ak, A%) < 
A. If packets are routed using the backpressure routing algo¬ 
rithm, then there exist finite constants T and R such that 
the following happens: 


1- For some t < T, Q„(t) > R for all n £ Ak, and 

2- For all t, Q„(t) < R for n £ A^- 

Proof. We will prove the two claims separately. To 
prove the first claim we will use the fact that the network 
is overloaded and bottlenecked at the cut {Ak,A%)- We 
will prove the second claim using the fact that the num¬ 
ber of packets that arrive into A^j. in each time-slot is upper- 
bounded by /™“, and any cut in the network has a capacity 
larger than or equal to The detailed proofs for both 

claims are given in the Appendix G. □ 

Each node n has a threshold-based smallest min-cut de¬ 
tection mechanism. When we start using a particular DAG 
Dk, in each time-slot, we check whether the queue crosses 
a prespecified threshold Rk- Any queue that crosses the 
threshold gets marked as overloaded. After using the DAG 
Dk for Tk timeslots, all the nodes that have their queue 
marked overloaded form the set Ak- When the time Tk 
and threshold Rk are large enough, the cut {Ak,Ak) is the 
smallest min-cut as proven in Theorem 3. After determin¬ 
ing the smallest min-cut, an individual node can perform a 
link reversal by comparing its queue’s overload status with 
its neighbor’s. All the links that go from a non-overloaded 
node to an overloaded node are reversed to obtain Dk+i- 
The complete LFBP algorithm is given in Algorithm 3. 


Algorithm 3 LFBP (Executed by node n) 

1: Input: sequences {Tk},{Rk}, unique ID n 
2: Generate initial DAG Dq by directing each link {n,j} 
to (n, j) if n < j, to (j,n) if j > n. 

3: Mark the queue Qn as not overloaded 
4: Initialize t 0, fc 0 
5: while true do 

6: Use BP to send/recive packets on all links of node n 

7: if {Qn{t) > Rk) then 

8: Mark Q„ as overloaded. 

9: end if 

10 : 

11 : 

12 : Tk^Tk-1 

13: if Tfe = 0 then 

14: Reverse all links {j, n) such that Qj is not over¬ 

loaded and Qn is overloaded. 

15: fc A: -f 1 

16: Mark Qn as not overloaded 

17: end if 

18: end while 


Corollary 1. Suppose the traffic demand is feasible in 
G, i-C-, A < Then, the LFBP algorithm (Algorithm 

3) will find a DAG, whose maximum flow supports A, in a 
finite number of iterations- 

Proof. Theorem 3 shows that LFBP identifies the small¬ 
est min-cut (A, A'’) for the DAG Dk- Lemma 6 shows that 
A is the set of overloaded nodes, and A'’ is the set of non- 
overloaded nodes in a flow allocation that induces the lex¬ 
icographically minimal overload vector. LFBP reverses the 
links going from A'’ to A, which is also the reversals per¬ 
formed by the link reversal algorithm (Algorithm 1). Hence, 
by Theorem 1, LFBP obtains the DAG that supports A. □ 






Good choices for the thresholds Tk and Rk are topology 
dependent. When the value of Rk is too small, nodes that 
are not overloaded might cross the threshold producing a 
false positive. If the value of Rk is large but Tk is small, 
the overloaded nodes might not have enough time to de¬ 
velop the backlog to cross Rk which produces false negatives. 
Hence, a good strategy is to choose a large Rk so that the 
non-overloaded nodes don’t (or rarely) cross this threshold, 
then chose a large Tk such that the overloaded nodes have 
enough time to build the backlog to cross Rk- Optimizing 
these thresholds requires further research. Note that our 
algorithm performance degrades graciously with false posi¬ 
tives/negatives. Even when it detects the smallest min-cut 
incorrectly, the actions of the algorithm preserve the acyclic 
structure. Thus, in the subsequent iterations the algorithm 
can improve the DAG again. 

4.1 Algorithm modification for topology changes 

In this section we consider networks with time-varying 
topologies, where several links of graph G may appear or 
disappear over time. Although the DAG that supports A 
depends on the topology of G, our proposed policy LFBP 
can adapt to the topology changes and efficiently track the 
optimal solution.Additionally, the loop free structure of a 
DAG is preserved under link removals. Thus, if some of the 
links in the network disappear, we may continue using LFBP 
on the new network. 

To handle the appearance of new links in the network 
smoothly, we will slightly extend LFBP to guarantee the 
loop free structure. For a DAG Dk, every node n stores a 
unique state x„{k) representing its position in the topolog¬ 
ical ordering of the DAG Dk- The states are maintained 
such that they are unique and all the links go from a node 
with the lower state to a node with the higher state. When 
a new link {i,j} appears we can set its direction to go from 
i to j if Xi{k) < Xj{k) and from j to i otherwise. Since 
this assignment of direction to the new link is in alignment 
with the existing links in the DAG, the loop-free property is 
preserved. 

The state for each node n can be initialized using the 
unique node ID during the initial DAG creation, i.e. a;n(0) = 
n. Then whenever a reversal is performed the state of node 
n can be updated as follows: 

- _ ( x-n{k — 1) — 2*’A, if n is overloaded, 

( x-n{k — l), otherwise. 

Here, A is some constant chosen such that A > maxijgjv a;i(0) — 
Xj{0)- Note that this assignment of state is consistent with 
the way the link directions are assigned by the link rever¬ 
sal algorithm. The states for the non-overloaded nodes are 
unchanged, so the links between these nodes are unaffected. 
Also, the states for all the overloaded nodes are decreased 
by the same amount 2*’A, so the direction of the links be¬ 
tween the overloaded nodes is also preserved. Furthermore, 
the quantity —2*’A is less than the lowest possible state be¬ 
fore the fcth iteration, so the overloaded nodes have a lower 
state than the non-overloaded nodes. Hence, the links be¬ 
tween the overloaded and non-overloaded nodes go from the 
overloaded nodes to the non-overloaded nodes. 

In this scheme, the states x„ decrease unboundedly as 
more reversals are preformed. In order to prevent this, after 
a certain number of reversals, we can rescale the states by 
dividing them by a large positive number. This decreases the 


value of the state while maintaining the topological ordering 
of the DAG. The number of reversals k can be reset to 0, 
and a new A can be chosen such that it is greater than the 
largest difference between the rescaled states. 


5. COMPLEXITY ANALYSIS 

To understand the number of iteration the link-reversal 
algorithm takes to obtain the optimal DAG, we analyze the 
time complexity of the algorithm. 


Theorem 4. Let G be a vector of the capacities of all 
the links in E, and let I be the set of indices 1, 2,..., |i?|. 
Define 5 > 0 to be the smallest positive difference between 
the capacity of any two cuts- Specifically, S is the solution 
of the following optimization problem 


min Ca 

A,BCI 

as A 


bGB 


subject to: Ca > Cb- 

aSA bSB 


The number of iterations taken by the link reversal algorithm 
before it stops is upper bounded by —] , where 

is the max-flow of the undirected network. 

Proof. After each iteration of the link-reversal algorithm, 
either the max-flow of the DAG increases, or the max-flow 
stays the same and the number of nodes in the source side 
of the smallest min-cut increases (see Lemma 8 in the Ap¬ 
pendix). We can bound the number of consecutive iterations 
such that there is no improvement in the max-flow. In par¬ 
ticular, every such iteration will add at least one node to the 
source set. So, it is impossible to have more than | A| —2 such 
iteration. Hence, every \N\ iterations we are guaranteed to 
have at least one increase in the max-flow. 

Max-flow is equal to the min-cut capacity, and min-cut 
capacity is defined as the sum of link capacities. Say, the 
max-flow of DAG Dk+i is greater than that of Dk- Let A be 
the set of indices (in the capacity vector G) of the links in 
the min-cut of Dk+i , and B be the set of indices of the links 
in the min-cut of Dk- This choice of A and B forms a feasible 
solution to the optimization problem given in the theorem 
statement. Since the optimal solution S lower bounds all the 
feasible solutions in the minimization problem, the increase 
in the max-flow must be greater than or equal to <5. 

Every \N\ iteration the max-flow increases at least by S. 
Hence, the DAG supporting the max-flow is formed 

within \\N\f’^’^^/5~\ iterations. □ 


Corollary 2. In a network where all the link capacities 
are rational with the least common denominator D G N, the 
number of iterations is upper bounded by {\N\Df’^’^^). 

Proof. Since the capacities are rational we can write the 
capacity of the i*’*' link as a = ^, where A/) is a natural 
number. From the definition of S in Theorem 4, we get 5 to 
be the value of the following optimization problem: 

min i ( X] -Yo - Aff, ] 
a,bciT> \ / 

\aeA b<^B / 

subjsct to. ^ ^ 2 , ^ Y^Mb. 

aSA beB 

All the A/(,) are integers, so to satisfy the constraint we must 
have the difference yA/a — A/i, > 1. Hence 5 > 



Using this value of S in Theorem 4, we can see that the 
number of iterations is upper bounded by {\N\'Df’^^^). □ 

Corollary 3. In a network with unit capacity links, the 
number of iterations the link-reversal algorithm takes to ob¬ 
tain the optimal DAG is upper bounded by |A'^||i?|. 

Proof. The max-flow < |i5|. So, by Corollary 2, 

the number of iterations is upper bounded by |A'^||i5|. □ 

We conjecture that these upper bounds are not tight, and 
finding a tighter bound will be pursued in the future re¬ 
search. We simulated the link reversal algorithm in 50,000 
different Erdos-Renyi networks (p — 0.5) of sizes 10 to 50 
with randomly assigned link capacities. The link reversal al¬ 
gorithm started with a random initial DAG. We found that 
it took less than 2 iterations on average to find the optimal 
DAG. 

A worst case lower bound for the number of iteration is 
lA^I. This lower bound can be achieved in a line network 
where the initial DAG has all of its links in the wrong direc¬ 
tion. 


6. SIMULATION RESULTS 

We compare the delay performance of the LFBP algorithm 
and the BP algorithm via simulations. We will see that the 
network with the LFBP routing has a smaller backlog on 
average under the same load. This shows that the LFBP 
algorithm has a better delay performance. We consider two 
types of networks for the simulations: a simple network with 
fixed topology, and a network with grid topology where the 
links appear and disappear randomly. 

6.1 Fixed topology 

We consider a network with the topology shown in Fig¬ 
ure 4(a). The edge labels represent the link capacities. 
The undirected network has the maximum throughput of 
15 packets per time slot. Figure 4(b) shows the initial DAG 
Do- Instead of running the initial DAG algorithm of Sec¬ 
tion 3.1, here we choose a zero throughput DAG to test 
the worst-case performance of LFBP. The arrivals to the 
network are Poisson with rate A = 15p , where we vary 
p = .5, .55,..., .95. For the LFBP algorithm, we set the over¬ 
load detection threshold to Rk = 60 for all n,k. To choose 
this parameter, we observed that the backlog buildup in nor¬ 
mal operation rarely raises above 60. We also choose the 
detection period T\ — 150 and T*, = 50 for all fc > 1. This 
provides enough time for buildup, which improve the accu¬ 
racy of the overload detection mechanism. 

We simulate both algorithms for one million slots, using 
the same arrival process sample path. Figures 4(c) - 4(e) 
show the various DAGs that are formed by the LFBP algo¬ 
rithm at iterations k = 1, 2,3. We can see that the nodes in 
the smallest min-cut get overloaded and the link reversals 
gradually improve the DAG until the throughput optimal 
DAG is reached. 

Figure 5 compares the total average backlog in the net¬ 
work for BP and LFBP, which is indicative of the average 
delay. A significant delay improvement is achieved by LFBP, 
for example at load 0.5 the average delay is reduced by 66% 
We observe that the gain in the delay performance is more 
pronounced when the load is low. In low load situations, the 
network doesn’t have enough “pressure” to drive the packets 
to the destination and so under BP the packets go in loops. 



chosen so that LFBP 
requires several itera¬ 
tion to reach the opti¬ 
mal. 



(e) The optimal DAG. 


Figure 4: Figure (a) depicts the original network. 
Figures (b)-(e) are the various stages of the DAG. 
The red nodes represent the overloaded nodes, and 
the dashed line shows the bonndary of the over¬ 
loaded and the non-overloaded nodes. 



Load 


Figure 5: Average backlog in the network (Fig. 4(a)) 
with fixed topology for the Loop Free Backpressnre 
(LFBP) and the Backpressnre (BP) algorithms. 


6.2 Randomly changing topology 

To understand the delay performance of the LFBP al¬ 
gorithm on networks with randomly changing topology, we 
consider a network where 16 nodes are arranged in a 4 x 4 
grid. All the links are taken to be of capacity six. For the 
LFBP algorithm, we choose a random initial DAG with zero 
throughput shown in Figure 6. The source is on the upper 
left corner (node 1) and the destination is on the bottom 
right (node 16). 
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Figure 6: Initial DAG for the LFBP algorithm cho¬ 
sen so that the LFBP needs several iterations to 
reach the optimal DAG. All the links have capacity 
six. 


In the beginning of the simulations all 24 network links 
are activated. At each time slot an active link fails with 
a probability 10“'* and an inactive link is activated with a 
probability 10“^. The maximum throughput of the undi¬ 
rected network without any link failures is 12. Clearly on 
average, each link is “on” a fraction ^ of the time, and 
thus the average maximum throughput of the undirected 
network with these link failure rates is x 12 = 10.9. The 
arrivals to the networks are Poisson with rate A = 10.9p, 
where p = .1, .2,..., .6. For the LFBP algorithm, the detec¬ 
tion threshold is set to Rk = 100 and the detection period is 
Tk = 30 for all n,k. These parameters were chosen so that 
there are several reversals before a topology change occurs 
in the undirected network. The simulation was carried out 
for a million slots. 

Figure 7 compares the average backlog of LFBP and BP. 
In the low load scenarios LFBP reduces delay significantly 
(by 85% for load = 0.1) even though the topology changes 
challenge the convergence of the link-reversal algorithm. As 
the load increases, both the algorithms begin to obtain a 
similar delay performance. 



Figure 7: Average backlog in the network with ran¬ 
dom link failures (Fig. 6) for the Loop Free Back¬ 
pressure algorithm and the Backpressure algorithm. 


7. MULTICOMMODITY SIMULATION 

We extend of the link reversal algorithm to the networks 
with multiple commodities. The multi-commodity algorithm 


is identical to the single commodity algorithm, with the ex¬ 
ception that we now use the multicommodity backpressure 
of [1]. Each node n maintains a queue Qn{t) for each com¬ 
modity y. Each commodity is assigned its own initial DAG. 
A pseudocode for the multicommodity LFBP that we used is 
given in Algorithm 4. An important direction for future re¬ 
search is to determine whether the claims proven for a single 
commodity in the previous sections extend to the multicom¬ 
modity case. 


Algorithm 4 Multicommodity LFBP (Executed by n) 

1: Input: sequences {Tk}, {Rk}, unique ID n 
2: For each commodity y, generate initial DAG Dq by di¬ 
recting {n,j} to (n, j) if n < j, to {j,n) if j > n. 

3: Mark all queues as not overloaded 
4: Initialize t 0, fc 0 
5: while true do 

6: Use Multicommodity BP to send/recive packets on 

all links of node n 

7; for all y do 

8: if > Rk) then 

9: Mark this Qn as overloaded. 

10: end if 

11: end for 

12 : 

13: 

14: Tk^Tk-l 

15: if Tfe = 0 then 

16: for all y do 

17: Reverse links {j,n) in if is not over¬ 

loaded and Qn is overloaded. 

18: end for 

19: k^kTl 

20: Mark all queues as not overloaded 

21: end if 

22: end while 


For the simulation, we consider a network arranged in a 
4x4 grid as shown in Figure 6. Each link has a capacity 
of 6 packets per time-slot. There are three commodities in 
the network defined by the source destination pairs (1,16), 
(4,13) and (5,8). For the LFBP algorithm, each commodity 
starts with the same initial DAG given in Figure 6. 

We use the arrival rate vector A™'*’' = [7.18,6.96,9.86], 
which is a max-flow vector for this network computed by 
solving a linear program. We scale this vector by various 
load factors p ranging from 0.1 to 0.9. The arrivals for each 
commodity i is Poisson with rate pA™“. In the beginning 
of the LFBP simulation, [500/pJ dummy packets are added 
to the source of each commodity. This is helpful in low load 
cases because it forces the algorithm to find a DAG with 
high throughput, and avoids stopping at a DAG that only 
supports the given (low) load. Rk was chosen to be 50 and 
Tk = 50 for all k > 0. The simulation was executed for 
500,000 time-steps. 

Figure 8 shows the average backlog in the network for dif¬ 
ferent loads under backpressure and multicommodity LFBP. 
We can see that the LFBP algorithm has a significantly im¬ 
proved delay performance compared to backpressure. 

8. CONCLUSION 

Backpressure routing and link reversal algorithms have 






























Load 


Figure 8: Average backlog in a multicommodity net¬ 
work with fixed topology for LFBP and BP algo¬ 
rithms. 

been separately proposed for mobile wireless networks ap¬ 
plications. In this paper we show that these two distributed 
schemes can be successfully combined to yield good through¬ 
put and delay performance. We develop the Loop-Free Back¬ 
pressure Algorithm which jointly routes packets in a con¬ 
strained DAG and reverses the links of the DAG to improve 
its throughput. We show that the algorithm ultimately re¬ 
sults in a DAG that yields the maximum throughput. Addi¬ 
tionally, by restricting the routing to this DAG we eliminate 
loops, thus reducing the average delay. Future investiga¬ 
tions involve optimization of the overload detection param¬ 
eters and studying the performance of the scheme on the 
networks with multiple commodities. 

9. REFERENCES 

[1] L. Tassiulas and A. Ephremides, “Stability properties 
of constrained queueing systems and scheduling for 
maximum throughput in multihop radio networks,” 
IEEE Transactions on Automatic Control, vol. 37, no. 
12, pp. 1936-1949, December 1992. 

[2] L. X. Bui, R. Srikant and A. Stolyar, “A novel 
architecture for reduction of delay and queueing 
structure complexity in the back-pressure algorithm,” 
lEEE/ACM Transactions on Networking, vol. 19, 
no. 6, pp. 1597-1609, December 2011. 

[3] M. J. Neely, E. Modiano and G. E. Rohrs, “Dynamic 
power allocation and routing for time varying wireless 
networks,” IEEE Journal on Selected Areas in 
Communications, Special Issue on Wireless Ad-hoc 
Networks, vol. 23, no. 1, pp. 89-103, January 2005. 

[4] E. Gafni and D. Bertsekas, “Distributed algorithms for 
generating loop-free routes in networks with 
frequently changing topology,” IEEE Transactions on 
Communications, vol. 29, no. 1, pp. 11-18, January 
1981. 

[5] V.D. Park and M.S. Corson, “A highly adaptive 
distributed routing algorithm for mobile wireless 
networks,” INFOCOM, 1997. 

[6] L. Georgiadis and L. Tassiulas, “Optimal overload 
response in sensor networks.” IEEE Transactions on 
Information Theory, vol.52, no. 6, pp. 2684-2696, June 
2006. 


[7] H. Xiong, R. Li, A. Eryilmaz and E. Ekici, 
“Delay-aware cross-layer design for network utility 
maximization in multi-hop networks.” IEEE Journal 
on Selected Areas in Communications, vol. 29, no. 5, 
pp. 951-959, May 2011. 

[8] L. Ying, S. Shakkottai, A. Reddy and S. Liu, “On 
combining shortest-path and backpressure routing over 
multihop wireless networks,” lEEE/ACM Transactions 
on Networking, vol. 19, no. 3, pp. 841-854, June 2011. 

[9] P.-K. Huang, X. Lin, and C.-C. Wang, “A 
low-complexity congestion control and scheduling 
algorithm for multihop wireless networks with 
order-optimal per-flow delay,” lEEE/ACM Trans, on 
Networking, vol. 21, no. 2, pp. 2588-2596, April 2013. 

[10] M. J. Neely, “Stochastic network optimization with 
application to communication and queueing systems,” 
Morgan & Claypool, 2010. 

[11] L. R. Ford and D. R. Fulkerson, “Maximal flow 
through a network,” Canadian Journal of 
Mathematics, 8: 399, 1956. 

[12] L. Georgiadis, P. Georgatsos, K. Floros, and 

S. Sartzetakis, “Lexicographically optimal balanced 
networks,” IEEE/ACM Transactions on Networking, 
vol. 10, no. 6, pp. 818-829, December 2002. 

[13] L. Huang and M. J. Neely, “Delay reduction via 
Lagrange multipliers in stochastic network 
optimization,” IEEE Transactions on Automatic 
Control, vol. 56, no. 4, pp. 842-857, April 2011. 

[14] L. Georgiadis, M. J. Neely and L. Tassiulas, “Resource 
allocation and cross-layer control in wireless 
networks,” Foundations and trends in networking. 

Now Publishers Inc, 2006. 

APPENDIX 
A. LEMMA 7 

Lemma 7. Consider a DAC Dk with source node s, des¬ 
tination node d, and arrival rate X. Let Ak be the set of 
overloaded nodes under the flow allocation (fij) that yields 
the lexicographically minimum overload vector. If \Ak\ > 0, 
then [Ak,A’/) is a min-cut of the DAG Dk. 

Proof Proof of Lemma 7. First we show that {Ak,A’/) 
is a cut, i.e., the source node s € Ak and the destination 
node d € A|. The destination node d has zero queue over¬ 
load rate qs = 0 because it does not buffer packets; hence 
d £ Ak. We show s £ Ak hy contradiction. Assume s ^ Ak. 
The property (8) shows that there is no flow going from Al 
to Ak, i.e., 

E = 0- 

(iJ)eEk:ieAl,jeAk 

The flow conservation equation applied to the collection Ak 
of nodes yields 

^ ' 9n = ^ ' fin — ^ ^ fnj 

neAk (i,n)eEk:i€A'^,neAk (n.j) 6B),: n6 , j £ A = 

E ^ 0’ 

(n,j)eEk: neAk,j€A<j_ 

which contradicts the assumption that the network is over¬ 
loaded (i.e., \Ak\ > 0). Note that in the above equation A 
does not appear because of the premise s ^ Ak. 





By the max-flow min-cut theorem, it remains to show that 
the capacity of the cut {Ak,A%) is equal to the maximum 
flow of the DAG Under the flow allocation {fij) 
that induces the lexicographically minimal overload vector, 
the throughput of the destination node d is the maximum 
flow (see Lemma 1). It follows that 


A - E 9* = 

A - E b-i 

(10) 

i^N 

iCAk 


E 

fij 

(11) 

(» iSAfc 



E 

Cij = capk{Ak,Al)- 

(12) 


46-4- 



where (10) uses (7) and qi = 0 for all nodes i ^ Ak, (11) fol¬ 
lows the flow conservation law over the node set Ak, and (12) 
uses the property (9) in Lemma 1. □ 

B. PROOF OF LEMMA 6 

Proof of Lemma 6. Lemma 7 shows that (AkjA'f.) is 
a min cut of the DAG Dk- It suffices to prove that if 
there exists another min-cut i.e., Ak ^ B and 

ca.Pk{Ak, Al) = ceiPk{B, B‘^), then Ak C B. The proof is 
by contradiction. Let us assume that there exists another 
min-cut {B, £?“) such that Ak (f. B. We have the source node 
s € Ak n B and the destination node d £ Alt'] B’^. Consider 
the partition {C,D,E,F} of the network nodes such that 
C ^ AkDB, D = Ak\B, E = B\Ak and F = N\{Ak U B) 
(see Fig. 9). Since Ak B and Ak yf B, we have \D\ > 0. 
Also, we have s £ C and d £ F. Let (fij) be a flow allo- 



Figure 9: A partition of the node set N where Ak = 

C U D and B = C L) E. 

cation that yields the lexicographically minimum overload 


vector in Dk- Properties (8) and (9) show that 

fij = Cij, Wi £ Ak, j £ Al, (13) 

fij =0, Vi € Ak, j £ Ak- (14) 

The capacity of the cut (B, B“) in the DAG Dk, defined 
in (1), satisfies 

caPfe(B,B“) = ca.Pk{B,D) + capkiB,F), (15) 
where B‘^ = DuF. Under the flow allocation {fij), we have 

caPfc(B,D)= Cij> Y /l- 


(i,j)GE}.:i£B,j£D ,j £D 

( 16 ) 


Applying the flow conservation equation to the collection of 
nodes in D yields 

E fij>Y^^+ E /d- (17) 

(-i,j)^Ek-i€B,j£D ieD {i,j)eEf,-.ieD,jeF 

In (17), the first term is the sum of incoming flows into 
the set D; notice that there is no incoming flow from F to 
D because of the flow property (14). The second term is 
the sum of queue overload rates in D. The last term is a 
partial sum of outgoing flows leaving the set D, not counting 
flows from D to B; hence the inequality (17). From the flow 
property (13), the outgoing flows from the set D to F satisfy 

E = E 

(i,j)eEk:ieD,jeF ^iJ)^Ek:ieD,jeF 

Combining (15)-(18) yields 
caPkiB, B‘') = ca.Pk{B, D) + capj,(B, F) 

^E^*+ E Cij -£ca.Pk{B,F) 

(i,j)eEk:ieD,j£F 

> Y Cij + cap k{B,F) 

(i,j)€Ek:ieD,j€F 

= caPkiAkU B,F), (19) 

where the second inequality follows that all nodes in D are 
overloaded and > 0 for all n £ D. Inequality (19) shows 
that there exists a cut {Ak UB,F) that has a smaller capac¬ 
ity, contradicting that {B, B‘^) is a min-cut in the DAG Dk- 
Finally, we note that the partition {Ak,A).) is unique be¬ 
cause the lexicographically minimal overload vector is unique 
by Lemma 1. □ 

C. PROOF OF THEOREM 3 

Proof of the first claim. First we will show that the 
queue at the source Qs{f) crosses any arbitrary threshold 
Ri- We know that for some node n £ Ak, Qn{t) —>■ cxa as 
t —>■ oo because the external arrival rate to the source s £ Ak 
is larger than the rate of departure from set Ak, i.e. A > 
cap{Ak,A%)- The backpressure algorithm sends packets on 
a link (i,j) only if Qi{t) > Qj{t)- Hence, at any time-slot if a 
node b yf s has a large backlog, then one of its parents p must 
also have a large backlog. Qp can be slightly smaller than Qk 
because Qb might also receive packets from other nodes at 
the same time-slot. Specifically, Qp{t) > Qb{t + 1) — JDi Cib- 
Performing the induction on the parent of p we can see that 
the source node must have a high backlog when any node in 
Ak develops a high backlog. Note that the network is a DAG 
and the node n received packets form the source to develop 
its backlog, so the induction much reach the source node. 
Hence, when Qb{Ti) ^ Ri, Qs{t) > Ri for some t <Ti- 
Now we will show that every node in Ak crosses the thresh¬ 
old R- Let Bi C Ak be the set of nodes such that Q„ {t) > Ri 
for some time t <T\- We showed that s £ B\- We will show 
that when B\ Ak, there exists some set Bj, such that (i) 
Bj C Bj, and (ii) for every node n £ Bj, Qn{t) > R 2 for 
some t KTj- Here, R 2 and T 2 are large thresholds. 

Assume Bi yf A*. Let Gi = Afe\Bi, i.e all nodes in Ci 
haven’t crossed the threshold R\ until time T\- Let cbiCi 
be the total capacity of the links going from B\ to Ci, and 
cCiAg be the total capacity of the links going from Ci to A%- 
We have cbiCi > cciA| because {Ak,A'j.) is the smallest 




min-cut (see Figure 10). When the backlogs of the nodes 
of Bi are much larger than the nodes of Ci, the nodes in 
Cl receive packets from Bi at the rate of cbiCi packets per 
time-slot, and no packets are sent in the reversed direction. 
The rate of packets leaving the nodes in C is upper bounded 
by CbiAI which is smaller than the incoming rate. Hence, at 
least one node n' £ C must collect a large backlog, say larger 
than 7?2 < Ri- So, each node in the set B 2 = B\ U{n'} have 
a backlog larger than R 2 at some finite time T 2 . 



Figure 10: Let {Ak,A%) be the smallest min-cut. We 
showed that s £ B\. Say, cciA‘= > cbiCi then the cut 
(Bi, Bf) has the capacity of CbiA<^+CbiCi < ccLp{Ak, A^). 
This contradicts the assumption that {Ak,Al) is the 
smallest min-cut. So, cciA^ < cbiCi- 

Now using induction we can see that for Bm where m < 
l^fcl, Bm = and all the nodes in Bm cross a threshold 
R = min{Bi,..., Rm} by time T — max{ri, ..., Tm}- □ 

Proof of the second claim. We will use the following 
fact to prove this claim: for any subset of nodes S, if the 
number of packets entering S is lower than or equal to the 
number of packets leaving S on every time-slot, then the 
total backlog in S doesn’t grow. So, the backlog in each 
node of S is bounded. 

Assume a node b develops a backlog Qb{t) > Ri- Here Ri 
is a chosen such that 

Ri = \Al\ Cij + max Q„{0). 

ijeA% 

Consider a subset B of Al such that for every node i £ B 
and j £ C = A%\B, {Qiit) — Qj{t)) > Cij. The sets B and C 
must be nonempty because Qb{t) is large and Qdit) is zero, 
that IS, h £ B and d £ C. Note that backpressure doesn’t 
send any data from C to B. 

Let Cab be the capacity of the links going from A to B, 
and let cbc be the capacity of the links going from B to C. 
So, the number of packets entering B at timeslot t is upper 
bounded by cab- The number of packets leaving B is equal 
to cbc- Since {A,A‘^) is the smallest min-cut, cab < cbc 
(see Figure 11). Hence, the number of packets entering B 
is less than or equal to the number of packets leaving it at 
time t. 

Therefore as soon as one of the nodes crosses threshold 
Bi, the sum backlog becomes bounded. We can choose a 
threshold B S> Bi such that this threshold is never crossed 
by any nodes in Al- 



Figure 11: Let {Ak,A%) be the smallest min-cut. We 
showed that d£ C. Say, cab > cbc then the cut (B U 
Ak, (B U AkY) has the capacity of cbc + ca^c < cab + 
CAf.c = cap{Ak, A%). This contradicts the assumption 
that {Ak,Al) is the smallest min-cut. So, cab < cbc- 

D. LEMMA 8 

Lemma 8. Consider the case when X > The link 

reversal algorithm is applied on DAG Dk to obtain Dt+i- 
Let {Ak,A%) and {Ak+i, A%^i) be the smallest min-cuts of 
Dk and Dk+i respectively. Then, either capk{Ak, Al) > 
capk+i{Ak+i,Al+i), or capk{Ak, A%) = capk+i{Ak+i, AIj^Y) 
and \Ak+i\ > \Ak\ 



Figure 12: Here L represents the sum of the capaci¬ 
ties of the links going from one partition to the next 
in the DAG Dk, and represents the sum of the 
link capacities in the DAG Dk-ei- For example, 19 
and 19' represent the links that go from {Ak U Afe+i)'’ 
to {Ak n Afc+i) in DAGs Dk and Dk+i respectively. 

Proof. Consider the partitioning of the nodes as shown 
in Figure 12. For i = 1,..., 12, h represents the sum of the 
capacities of the links going from one partition to the next in 
the DAG Dk, and Z' represents the sum of the link capacities 
in the DAG Dk+i- The capacities of the smallest min-cut, 
before and after the reversal are given by 

capk{Ak,Al) = I 2 + I 5 + ho + h 2 and 

capk+i{Ak+i, Afc_|_i) = I 4 , h + ho + ^11 

respectively. Note that only the links that are coming into 
Ak are different in Dk and Dk+i- So 

li = I'i for i = 3,4,7,8,10,12. 


□ 


( 20 ) 






Because of the reversal there are no links coming into Ak in 
the DAG Dk+i’- 

( 21 ) 

After the reversal, the incoming links to Ak become outgoing 
from Ak, 

lio ~ ^10 + ^ 9 - ( 22 ) 

(Corresponding equations for I 2,15 and l'i 2 are omitted be¬ 
cause they are not necessary for the proof). Since {Ak,A%) 
is a min-cut, 

h < h- (23) 

This is true because otherwise the cut [Ak U Ak+i, (Ak U 
Ak+i)‘^) in the DAG Dk has a smaller capacity then the min 
cut {Ak,AkY. Specifically, let us assume I 5 > I 7 . Then, we 
get the contradiction: 

capk{Ak U Ak+i, (Ak U Ak+i)‘^) = 12 + 1? + ho 

< h + I 5 + ho + h 2 
= capk{Ak,AkY 

First we will show that if Ak\Ak+i Y 4‘^ then the capacity 
of the DAG must have increased. The proof is by contradic¬ 
tion. 

Let us assume that the throughput didn’t increase. So, 
capk{Ak,Al) > capk+i{Ak+i, Al+Y) 


II 

-b 

-b 

0 

-b 


= /4 + ^7 + lio + 0 

(24) 

> /4 + ^5 + lio 

(25) 

= capk{Ak n Ak+i, {Ak n Ak+iY). 

(26) 


(24) is follows from (20) and (21), and (25) follows from (23). 
Since Ak\Ak+i 7 ^ 0 by assumption, \Ak\ > \Ak n Afc+i|. 
This leads to a contradiction, because in DAG Dk the cut 
[Ak n Ak+i, [Ak n Ak+iY) is smaller than the smallest min- 
cut [Ak,Al). Hence, capk[Ak, Al) < capk+\[Ak+i, AIj^Y■ 
Next, we will consider the case Ak\Ak+i = <j>- Using (23), 

capk[Ak, Ak) = + Zio < h + ho- 

In this situation, we again have two cases. First, if Ak = 
Ak+i we know that Z(o > ho and Z7 = 0. Hence, capk[Ak, Al) < 
ho = co,Pk+i[Ak+i, Ak+Y)- 

Second, if Ak C Ak+i, then \Ak\ > |Afc+i| and 

ho > ho- (27) 

Using (20) and (27) capk[Ak,Al) < capk+i[Ak+i, Al^Y- □ 



